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A non-supersymmetric grand unified theory can exhibit a "radiative fermion mass hierarchy", 
in which the heavier quarks and leptons get mass at tree level and the lighter ones get mass from 
loop diagrams. Recently the first predictive model of this type was proposed. Here it is analyzed 
numerically and it is shown to give an excellent fit to the quark and lepton masses and mixings, 
including the CP phase violating phase 5c km- A relation between the neutrino angle 613 and the 
atmospheric neutrino angle is obtained. 
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, The masses of the known quarks and leptons exhibit a large hierarchy. This has suggested to many theorists 
' [l| that the light fermion mass hierarchy could be "radiative", i.e. that the lightest fermions get mass from loop 
' diagrams, while the heaviest get mass at tree level. In the early 1980's several papers showed that such an idea 
can be implemented naturally in the context of non-supersymmetric grand imified theories (GUTs) 0|. In models of 
this type, the radiative masses come from loop diagrams containing virtual GUT-scale particles. That is why such 
■ models must be non-supersymmetric: otherwise, the loops would be suppressed by 0{Mg^gY/^GUT) "^^^ *° 
Oh' non-renormalization theorems of supersymmetry. 

^ ^ In a recent paper Q a very simple non-supersymmetric 5*0(10) model with a radiative hierarchy was proposed. 
^ One thing that allows this model to be so simple is precisely that its hierarchy is radiative. The point is that terms 
, have to exist in the lagrangian corresponding to the larger elements of the quark and lepton mass matrices, but not 
to the smallest elements, since they arise automatically from loops. The simplification can be seen by comparing the 
model of [J] to the supersymmetric 50(10) model on which it was based, which had a non-radiative hierarchy 
^ . That earlier model had a somewhat larger particle content and more Yukawa terms. 
fSJ ' Models with radiative hierarchies are also simpler in another way: in them it is not necessary to introduce ad 
T-H , hoc very small dimensionless parameters to account for the fermion mass hierarchies, since they are automatically 
' accounted for by the loop factors I/IGtt^. Despite radiative hierarchy models being able to have a simpler structure, 
, one might think they would be less predictive, since loop diagrams tend to depend on many parameters. However, 
[ — the model proposed in [^l shows that this need not be the case. In that paper it was shown that the model gives a 

, qualitatively realistic pattern of quark and lepton masses and mixings with only 9 parameters. 
00 . While the non-supersymmetric model proposed in Q is economical and qualitatively realistic, the analysis in that 
paper was not sufficient to establish that it is realistic quantitatively. In particular, several issues were not addressed. 
First, it was not specified what the sequence and scales of breaking were of 50(10) down to the Standard Model group 
GsM (which must, of course be consistent with proton decay bounds and unification of gauge couplings). Unless that 
^ , is done, the renormalization-group running of the quark and lepton masses needed for a global fit of parameters cannot 
' be performed. Second, the forms of the mass matrices given Q were derived under the assumption that certain SU{5)- 
breaking effects could be ignored. However, as will be seen, this assumption is not necessarily consistent with the 
pattern of 50(10) breaking that needs to be assumed in order to satisfy the constraints of gauge coupling unification 
and proton decay. 5?7(5)-breaking effects will turn out to modify significantly the forms of the quark and lepton mass 
matrices given in [3|. Third, the mb/rrir ratio is problematic in the version of the model discussed in [5]. In that 
model, because certain 5J7(5)-breaking effects were treated as negligible, the classic prediction to" = m° was obtained. 
(A superscript '0' indicates throughout this paper quantities evaluated at the GUT scale.) This is well-known to give 
a fairly good fit in supersymmetric models for certain values of tan /3 Q ; but in non-supersymmetric models it results 
in a prediction of rnb/rrir at low energies that is typically too large by at least 30% 6]. Fourth, the ratio rris/mb 
is predicted in the version of the model given in to have the Georgi-Jarlskog value ^m^/mT at the GUT scale 
However, lattice calculations Q have suggested that to^ is significantly smaller than previous estimates of it, and the 
best fit value is now somewhat smaller than the Georgi-Jarlskog prediction. 

In this paper, we address all these issues. The paper is organized as follows. In section 2, the 5*0(10) model of Q 
is reviewed, and it is explained how both the tree-level and radiative contributions to the mass matrices arise, and 
why the resulting forms give a good qualitative description of the pattern of quark and lepton masses and mixings. In 
section 3, a breaking of 50(10) down to the Standard Model consistent with gauge coupling unification and proton 
decay bounds is specified. In section 4, the effect of this pattern of symmetry breaking on the quark and lepton mass 
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matrices is discussed and it is shown that forms somewhat different from those given in [3| result. In section 5, the 
results of a global numerical fit to the quark and lepton masses and mixings is given. An excellent fit is found to the 
quark and lepton masses and mixings, including the CP phase 5c km- A relation between the neutrino angle Oiz and 
the atmospheric neutrino angle is obtained. 



II. THE MODEL 



The model proposed in [3|, whose predictions we analyze in detail in this paper, is a non-supersymmetric with 
unified group 5*0(10). In it the tree-level mass matrices of quarks and charged leptons are generated by only three 
effective Yukawa operators 



01 = I63I63IOH 

02 = 162l63l0H45H/AfGC/T (1) 

03 = (ca6a6.H)(163l6'^)/MGt/T , i = 1, 2 

In O3, the factors in parentheses are contracted into lO's of S'O(IO). The loop-level elements in the mass matrices 
arise very simply from the tree-level elements, as will be seen later. The three operators given in Eq. (I) do the 
following things: Oi gives the 33 elements of the mass matrices, i.e. the masses of the third family. O2 and O3 
generate the masses of the second family and its mixing with the third family (i.e. Vcb and Oaun), and also 6 sol- The 
masses of the first family and its mixings come from loops. 

The operators in Eq. (1) come from integrating out some "extra" vectorlike fermion multiplets, consisting of an 
SO{lQ) vector and a spinor-antispinor pair. Thus the complete fermionic content of the model comprises the following 
(left-handed) multiplets: 16j=i,2,3 -I- (16 -I- 16 -(- 10). The Dirac mass matrices of the up-type quarks, down-type 
quarks, charged leptons, and neutrinos (denoted by Mjj, M]j, Mi^, and M^, respectively) arise from the following set 
of Yukawa terms in the lagrangian: 

jCvuk = Afi6(T6 16)-f Mio(10 10). 

+ a(T6 163)45ff + E.=i2C*(10 16,)16,ff 

(2) 

+ /i33(163l63)10ff + /i2(16 162)10jf + hsiW 163)16^ 



-I- h{16 16)10'^:^. 

It is shown in ^ that this form of the Yukawa interactions is the most general allowed by a certain simple U{1) flavor 
symmetry, which will be denoted by U{1)f- The terms on the first line of Eq. (2) are the 0{Mgut) masses of the 
extra fermion multiplets; the terms on the second line contribute 0{Mgut) masses that mix those extra fermions with 
the three chiral families 16i; the terms on the third line generate the weak-scale SU{2)l x C/(l)F-breaking masses; 
and the last term is needed to give radiative masses to the first family. Higgs multiplets are denoted by the subscript 
H. The Higgs fields 16iH obtain vacuum expectation values (VEV) in the 1(16) direction. (The expression p(q) 
stands for a p multiplet of SU{5) contained in a q multiplet of 5*0(10).) The adjoint Higgs field 45h is assumed to 
obtain a VEV that is proportional to the 5*0(10) generator B — L (i.e. baryon number minus lepton number). 

The electroweak gauge symmetry SU{2)]^ (g) U{1)y is spontaneously broken by the Higgs multiplets denoted IOh, 
IO'h, and 16^ in Eq. (1), and, more specifically, by the neutral components of the Y/2 — —1/2 doublets contained 
in 5{10h), 5(10'^), and 5(16^), and the neutral components of the Y/2 = +1/2 doublets contained in 5(10^) and 
5(10'^). Of course, in the low-energy effective theory, which is just the Standard Model, there is only one Higgs 
doublet, which is some linear combination of these doublets (and their hermitian conjugates). 

According to [3], the mass matrices that result from the terms in Eq. (2) have the form 
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where mu = /i33(5(10//)) and niD = /i33(5(10^f )). (It will be seen in section 4 that GUT-symmetry-breaking effects 
modify these forms somewhat.) The convention here is that the mass matrices are multiplied from the left by the 
left-handed fermions and from the right by the right-handed fermions. 

The 33 elements of the mass matrices in Eq. (3) come simply from the term /i33(163 163)10//, as is usually the 
case in 5*0(10) models [9]. (This is just the operator Oi in Eq. (1).) 

The contributions to the 23 and 32 elements denoted by e come from integrating out the family-antifamily pair 
16 + 16. The antifamily 16 appears in two mass terms from Eq. (2), which can be combined as follows: 16(Mi6l6 -I- 
a(45/f)163). These terms have the effect of mixing the 16 with the I63. One linear combination of 16 and I63 
obtains an 0{Mgijt) mass, while the orthogonal combination (denoted by the index 3') remains light. (From now 
on, primed indices will be used to denote the light families that remain after the superheavy fermions have been 
integrated out.) Thus, the 16 with no index has some of the third light family mixed in with it; and the amount of 
this mixing is proportional to the VEV (45/f). As a result, the term ^-2(16 162)10/f from Eq. (2) leads to an effective 
operator of the form (163/162' )10i/45_f/A^G;7T, which is just the operator O2 of Eq. (1), which in turn produces the 
contributions denoted in Eq. (3) by e. Since (45//) oc B — L, the e contributions are 1/3 times as large for the quarks 
as for the leptons. 

The elements denoted by Ci and C2 arise in a similar fashion by integrating out the 5O(10)-vector multiplet of 
quarks and leptons, 10. This multiplet contains a 5 + 5 of SU{5). The 5(10) appears in several mass terms from Eq. 
(2), which can be combined as 5(10)[Afio5(10) + 2 Ci(l(16i//))5(16i)]. These terms have the effect of mixing 

the 5(10) with the 5(16i) and 5(162). One linear combination of these 5's obtains an 0{Mgut) mass, while the two 
orthogonal linear combinations are in the light families and denoted 5i' and 52'- Consequently, the 5(10) has mixed 
in with it some of 5i' and 52'. That means that the term /i3(10 163)16^ in Eq. (2) leads to effective mass terms 
of the form {Ciby + C252')103'm/). This is just the operator O3 of Eq. (1) and gives the terms denoted by the Ci 
in Eq. (3). These contributions appear only in Ml and Md, because 5's of SU{5) contain only charged leptons and 
down-type quarks. In both Q and [1] the Miq was assumed to be an explicit (and therefore S'{/(5)-invariant) mass, 
and therefore the same Ci appear in both Ml and M/j. 

At this point it should be noted that the expressions for the quark and lepton mass matrices given in Eq. (3) are 
approximate. The exact expressions involve factors, such as 1/a/1 -I- (a(45//) /Mig)^ and + Ci(16i//)/Afio)^, 

which are essentially just the cosines of angles describing the mixing between the extra fermions 16-1-16-1- 10 and the 
three chiral families 16^. If these mixing angles are small, their cosines are very close to one, and the mass matrices 
become insensitive to their values. This is an assumption that we make here (as in [^), as it reduces the number of 
parameters. However, there is no a priori reason to assume that these angles are extremely small. (Indeed, if they 
vanished, so would e.) If one of these angles were of order 0.25 radians, say, it would give 3% corrections to some of 
the elements of the mass matrices. 

The elements denoted by Sgi and Sh in Eq. (3) are necessary to make the mass matrices Ml and Md be of rank 3 
rather than rank 2, and so generate masses and mixings for the first family. As will be seen, in order to fit the first 
family masses and mixings these S are must be of order 10^^, whereas the other parameters appearing inside in the 
mass matrices in Eq. (3) turn out to be of order 1 (or, in the case of e, about 0.19). In 4], additional vectorlike quark 
and lepton fields besides those in Eq. (2) had to be introduced in order to generate these small (5's. In 3], however, 
it was noted that that the terms in Eq. (2) are enough to generate the S terms automatically by one-loop diagrams 
and also to explain why they are of order 10~^. 

The Sgi terms are given by the one-gauge-boson- loop diagram shown in Fig. 1(a). The gauge boson in this diagram 
is in a 10 of SU{5) (of course, it is in the adjoint 45 of 50(10)), so that it turns lO's of 5[/(5) into 5's and vice 
versa. That means that the small Sgi elements that couple 10^ to 53 (namely (M/j)i3 and {ML)3i) come from the 
large Ci elements that couple 5i to IO3 (namely (Mi)i3 and (M/))3i) So Sgi cx Ci. These diagrams were evaluated in 
3 neglecting certain 5C/(5)-breaking effects, giving the result that the same Sgi appear in Ml and Md, as given in 
Eq. (3). 
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Figure 1. The diagram in (a) sliows how a trcc-lcvcl mass for lOsSj (shown as a blob in the center) leads to a 
one-loop mass for 10i53: i.e. the 6gi elements arise radiatively from the Ci elements. The 10(45,) in the loop is a 
superheavy gauge boson. The diagram in (b) is more detailed and shows why the loop is finite. 



The diagram in Fig 1(a) superficially looks divergent. However, the accidental symmetry U{1)f of the terms in 
Eq. (1) makes {Md)iz and {Ml)zi vanish at tree level and guarantees that the loop is finite, as an exact calculation 
indeed shows. The finiteness of this diagram is more obvious if we write it in the form shown in Fig. 1(b). The 
calculation of these loops will be discussed in section 4. 

The 22 elements of Ml and Md (denoted 5h) arise from the one-Higgs-boson-loop diagram shown in Fig. 2. 
Whereas the one-gauge-boson-loop shown in Fig. 1 must exist if the tree-level masses in Eq. (3) exist, the diagram 
in Fig. 2 only exists if certain couplings not needed for the tree-level masses are present: namely, the last term in 
Eq. (2) (ft,(16 16) 10^) and a Higgs-mass term of the form lOjf lO/j. A diagram related by 5*0(10) to the one in 
Fig. 2 gives a 22 element for the up-quark mass matrix Mjj. However, if one supposes the contributions to {Mij)22 
and (M£))22 from these diagrams to be roughly comparable, then {Mu)22/{Mu)ss would be of order 10~^ and thus 
at most a few percent correction to rric. 
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Figure 2. A diagram showing how the 22 elements of the mass matrices can arise radiatively through Higgs-boson 
loops. 



Before getting into a more detailed discussion of the model, it is useful to explain how the structure of the matrices 
given in Eq. (3) explains qualitatively many of the features of the observed pattern of masses and mixings of the 
quarks and leptons. 

First, neglecting the S parameters (which are of order 10~^ because they come from one-loop diagrams) and the 
parameter e (which, though a tree- level effect, is somewhat smaller than 1), one has that all the mass matrices in Eq. 
(3) are of rank 1. In this approximation, = = a/1 -I- |Ci p -|- IC2P, where the superscript '0' denotes quantities 
evaluated at the GUT scale. The relation mJJ = is known to fit fairly well in supersymmetric grand unified models 
with certain values of tan (3. It works less well in non-supersymmetric grand unified models; however, this relation 
will be substatially modified when realistic symmetry breaking is taken into account in section 4. 

The large (i.e. 0(1)) off-diagonal elements Ci produce large mixing angles in the left-handed lepton sector and 
the right-handed quark sector. This is because they result from mixings of 5's of SU{5), which contain, of course, 
left-handed charged leptons and right-handed down-type quarks. Consequently, these elements produce large MNS 
neutrino mixing angles, but they do not produce large CKM mixing, since CKM mixing is of the left-handed not 
right-handed quarks. This is one of the basic ideas of so-called "lopsided" models [lol |. 

Moreover, the present model is "doubly lopsided" in the sense that both Ci and C2 are large [ll|. (In singly lopsided 
models C2 is large but not Ci.) This doubly lopsided structure can give a very simple explanation of the "bi-large" 
pattern of neutrino mixing angles in the following way. Diagonalizing the charged-lepton mass matrix Ml requires 
0(1) rotations to eliminate the elements Ci and C2. In particular, it requires a rotation by 9s = ta,n~^ {Ci / C2) in the 
1-2 plane followed by a rotation by 0a = tan~^ -\/|Cip + IC2P in the 2-3 plane. If, as here, the neutrino mass matrix 
is hierarchical, then the rotations required to diagonalize the charged lepton mass matrix Al^ give the dominant 
contributions to the MNS matrix, which is therefore approximately of the form 



(Cs Ss \ 

-CaSs CaCs Sa , (4) 
Sa^s SaCs Ca j 

where = sin^a, Cq = cos^aj = sinfl^, and = cos^^- This is the bi-large form, with ds being the solar neutrino 
angle and 9a being the atmospheric angle. When the effects of the small parameters e and the (5's are taken into 
account, there will be small shifts in VmnSi including a small non-zero ^13, which will be calculated in section 5. 

If one considers the effects of e, but still neglects the the i5's, one sees by inspection of the mass matrices in Eq. (3) 
that rric/mt is of order e^, whereas ms/mh, m^/rriT and Vet are all of order e. This corresponds to the experimental 
fact that rric/rnt — 0.0025, whereas ms/mi,, m^/mr: and Vcb are given respectively by 0.02, 0.06, and 0.04. In fact, it 
is easy to show from Eq. (2) that ~ | sin^ 9a and {ms/mb)^ ~ | sin^^ cos^q, so that Vcb ~ ms/rrib, in agreement 
with observation. One also easily sees that {nis/mb)^ ~ ^^ i+c^ (m^/mr)" — e ^^^^ , where C = + 
so that the ratio {ms/mb)'^ / {m^/mrf is approximately 1/3, the Georgi-Jarlskog prediction 0]. 
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The u quark is left massless by the matrices in Eq. (3). However, that is not a bad thing. Experimentally, mu/mt 
is of order 10~^, which is far smaller than the corresponding ratios rrid/mi, ~ 10"'^ and m.f,/mr ~ 0.3 x 10^"^. Thus, 
if TOrf and me arise at one-loop level, one would expect m„ to vanish at one-loop level. 

Before analyzing the structure of the quark and lepton mass matrices further, it is necessary to deal with the 
question of the pattern of breaking of 50(10) down to the Standard Model group Gsm — SU{3)c <8i SU{2) l®U 

III. THE BREAKING OF 5'O(10) 

If 5*0(10) broke at a single scale all the way to the Standard Model group Gsm, it would give the same prediction 
for the low energy gauge couplings as non-supersymmetric SU{5), which are known to be unsatisfactory. Moreover, 
as in non-supersymmetric SU{5), the proton lifetime would be too short. However, it is possible to get satisfactory 
gauge coupling unification and proton lifetime by assuming a two-stage breaking with the Pati-Salam group [l^] as 
the intermediate symmetry: 

50(10) ^ 5C/(4)c SU{2)l ® SU{2)r SU{3)c SU{2)l ® U{1)y 

The breaking of 50(10) to the Pati-Salam group at the higher scale Mq can be done by a 210//. The breaking of 
the Pati-Salam group at Mps is done by the VEVs of the adjoint and spinor Higgs fields, 45h and 16iH, which also 
contribute to the superheavy quark and lepton masses through couplings that appear in Eq. (2). 

In running the gauge couplings between Mq and Mps, the following matter multiplets contribute to the beta 
functions: (1) The quark and lepton multiplets, 16i, I62, I63, 16, 16, and 10. Since the masses of these multiplets 
are produced by coupling to the adjoint and spinor Higgs fields and not the 210/f , their splittings are of order Alps 
and can be ignored in running between Mq and Mps- (2) The Higgs multiplets IQih, I62//, 16^, 10/f, 10^, and 
three adjoint Higgs multiplets. For the Higgs multiplets too, except for the adjoints, it is assumed that the splittings 
are of order Mps and can be neglected in running between Mq and Mps- For the adjoints, however, we assume 
splittings of order Mq. Under the Pati-Salam group a 45jy decomposes to (15, 1, 1)h + (6, 2, 2)/f-f (1, 3, 1)h + (1, 1, 3)//. 
We assume that the color-singlet pieces of the adjoints get mass of order Mq and the color-non-singlet pieces get 
mass of order Mps- This is not unreasonable, since the renormalizable couplings of the adjoints to the 210/f produce 
splittings of order Mq between the color-singlet and color-non-singlct pieces. (One such term is (^H^^'^^^^)H^^'^^H^^^\ 
where the fundamental indices I,J,K,L take SU{2)l ® SU{2)p values.) Of course, the whole Higgs potential must 
be tuned to give the hierarchy between Mq and Mps, so different patterns of splittings are possible. It is in order to 
get a value of Mq large enough to be consistent with proton-decay limits, that we assume there are three split adjoint 
Higgs multiplets. More such adjoints would push Mq higher. Below the scale Mps, we assume just the Standard 
Model field content with one Higgs doublet. 

The results of the running are shown in Fig. 3. 




Figure 3. The running of the gauge couplings. The Pati-Salam scale is indicated by the vertical line. 



In the running, the input values used are \Mz) = 58.97, a-\Mz) = 29.61, &nda-\Mz) = 8.47 [H. The result of 
the running is that Mps = 4.79x 10^^ GeV. At Mps the Standard Model couplings have the values a^^{Mps) = 41.35, 



7 



a^^iMps) = 43.21, and a^^{Mps) = 38.55. The unification scale comes out to be Mq = 1.17 x lO^*^ GeV. and 
a^^{MQ) — 35.65. These values are consistent with present bounds on proton decay. (In this model, the dominant 
contribution to proton decay comes from dimension-6 operators produced by the exchange of gauge bosons of mass 
Ma- The Pati-Salam gauge bosons do not give proton decay.) The values of the gauge couplings plotted in Fig. 3 
are used in the running of the quark and lepton masses and mixing angles in section 5. 

One consequence of the fact that 50(10) is broken down to the Pati-Salam group at a high scale is that it 
makes more natural the assumption being made in this model that (45 jy) (x B — L. In the original supcrsymmctric 
version of the model [4] this assumption was justified by the fact that a 45h whose VEV is proportional to B ~ L 
is needed to implement to Dimopoulos-Wilczek mechanism (or "missing VEV mechanism") [l3i] for doublet-triplet 
splitting. However, in a non-supersymmetric SO{10) model, that mechanism does not work — and, in fact, doublet- 
triplet splitting must be achieved through fine-tuning [Tsj . The justification for the assumption made in Q that 
(45jj) (X B — L is therefore less clear. However, if 50(10) breaks to the Pati-Salam group at the high scale Mq in 
such a way that only the (1, 3, 1)h + (1, 1, 3)^ in the adjoints have the large mass Mq, as assumed, then the residual 
(15, 1, 1)h naturally obtains a VEV in the B — L direction (and, in fact can do no other without breaking the Standard 
Model group at a superheavy scale). 

IV. MODIFICATIONS TO THE MASS MATRICES 

In this section the implications for the mass matrices of the breaking of 50(10) down to the Pati-Salam group will 
be analyzed. 

In Eq. (3) the term Afio(10 10), as written, involves an explicit mass. This was the assumption made in Q. It 
is also possible, however, and just as simple, to suppose that this mass arises from the VEV of some Higgs fields 
that breaks 50(10). For example the term may come from (10 10)(45//)(45i/)/Mps or from (10 10)(54//). If Mw 
reflects the breaking of SU (5) then it is a matrix that gives different values when acting on the down quarks and 
on the charged leptons in the 50(10) 10 of fermions. Call its value for the leptons M and for the down quarks M'. 
One result of this splitting is that the e ntries Cj are no longer t he same in the mass matrices and Md . If one 
assumes, as before, that the quantities a/I -|- (SiCi(16i^f)/7\fio)^ are well approximated by 1, then one can write the 
Ci contributions to the mass matrices as Ci for Ml and fCi for AId, where / — M/M'. This 5O(10)-breaking effect 
is, as it were, optional. However, it is quite important for fitting m'^/m'^, which otherwise would be predicted to be 
1. The effects that will now be described are unavoidable consequences of the breaking of 50(10). 

5O(10)-breaking and 5[/(5)-breaking effects come into the loop contributions Sgi and Sh in several ways. Consider 
Sh first. If one examines the diagram in Fig. 2 closely, one finds that its contribution to Mu involves loops with 
scalars that can be either color triplets or color singlets, whereas its contributions to Ml involve only color-triplet 
scalars. When 50(10) breaks to the Pati-Salam group, the degeneracy between these two types of scalars is badly 
broken, which means that one cannot assume that Fig. 2 gives the same contribution to the two matrices. We 
therefore introduce a factor fn into the 22 element of Ml to reflect this fact. 

The case of the elements Sgi requires a more involved discussion. First, it must be noted that the vacuum expectation 
values (l(16iH)) and (45//) break the Pati-Salam group, and therefore must be no larger than Mps- Moreover, the 
masses Miq and Miq cannot be too much larger than these VEVs, since otherwise the entries Ci and e would be too 
small. Consequently, one can assume that all the superheavy fermion masses are much lighter than the scale Mq. 
This has implications for the loop diagrams in Fig. 1. Some of those diagrams contain gauge bosons whose mass is 
of order Mps and others contain gauge bosons whose mass is of order Mq- Because the fermions in those loops are 
much lighter than Mq, as just argued, the loops with 0{Mc) gauge bosons are suppressed relative to the loops with 
0{Mps) gauge bosons by a large factor and are therefore negligible. To see what this implies, one must look in more 
detail at the diagrams in Fig. 1 to see how they contribute to Ml and Mpi. 

In Fig. 4(a), is shown the contribution to Ml- In this diagram there are three possible values of the color index 
a, (or, equivalently, of the pair of color indices be on the superheavy gauge boson X*"^). For any of these values, the 
gauge boson converts a left-handed charged lepton into a left-handed down-type quark at one vertex, and a left-handed 
charged anti-lepton into a left-handed down-type anti-quark at the other vertex. That shows that the gauge boson is 
one of those of the Pati-Salam group SU (4)c, which make such transitions. (The Pati-Salam multiplet (4, 2, 1) unifies 
left-handed leptons and quarks, while the multiplet (4, 1,2) unifies left-handed anti-leptons and anti-quarks.) For all 
three values of a, therefore, the loop in Fig. 3(a) contains only gauge bosons whose mass is of order Mps- 

In Figs. 4(b) and 4(c) are shown the diagrams that contribute to Mp). In Fig. 4(b), the external quark has a 
color index of fixed value a, which determines uniquely the values of the color indices on the virtual gauge boson 
X'"^. This gauge boson is (as in Fig. 4(a)) one of those of the Pati-Salam group 5?7(4)c, as can be seen from the fact 
that it converts left-handed down quarks into left-handed charged leptons. In Fig. 4(c) there are two choices for the 
color index c on the gauge boson X^'^, since it is only required to be different from a. This gauge boson, however, is 
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obviously not one of those of the Pati-Salam group, since it converts a left-handed quark into a left-handed anti-quark 
(and a left-handed lepton into a left-handed anti-lepton), which are not unified together in the multiplets of the 
Pati-Salam group. Thus the gauge boson in Fig. 4(c) has mass of order Mq ^ Mps- The diagram in Fig. 4(c) is 
thus highly suppressed. From these considerations, one sees that the contribution to Ml has a factor 3 relative to the 
contribution to Mjj because of the color degeneracy in the loop. 



X 



be 



(3) 



V'(3)l 
= (■. 



(3) 



Fig. 4(a) 



X 



be 



= rf(i) 



-*(3) 



V'(3)a 



(3) 



Fig. 4(b) 



X 



Ic 



i'(i)b 



^(3) 



(3) 



^(3)a 

= d. 



(3) 



Fig. 4(c) 

Fig. 4 Gauge-boson-loop contribution to mass matrices. Subscripts in parentheses are family labels, a, 6, c are 
SU{i)c indices. 1,2 are SU{2)l indices. 



The gauge-boson-loop integrals can be written in a simple form if one makes the same approximation as before, 
namely ^Jl+'{^ifiJ^i^jJ]Y[MlQ^ = 1. In that case the gauge-boson- loop contributions to Mi, and are given by 



Here x = (Mg/M)"^ and x' = (Mg/M')'^ , where Mg is the mass of the colored "Pati-Salam" gauge bosons in Fig. 3(a) 
and 3(b). Recalling that / = M/M', one can rewrite these as 



9 



where F{x) = x^^^\nx/{x — 1). It happens that the function F(x) is very slowly varying for arguments of order 1. 
For example, F{1 + y) — I — -^y^ + + and F{2) = F{^) = 0.98. Thus to a very good approximation one may 
write 



{Mn)f/ = as. 



The mass matrices that result are 



(7) 



M 



N 






-e I mu, Mr 
e 1 






fCl fC2 



CiS 
Sh 1 + 6*2(5 I niD, 





fH&H 



(8) 



3Ci(5 e + 3C2(5 1 



V. FITTING THE FERMION MASSES AND MIXINGS 



The forms of the mass matrices in Eq. (8) are those given by the model at the scale Alps, since that is the mass of 
the superheavy fields that are integrated out to give these matrices. Below Mps, the model reduces to the Standard 
Model. One can therefore use the renormalization group equations (RGEs) of the one-Higgs-doublet Standard Model 
to run the measured quark and lepton masses and mixings from low scales up to Mps and then fit the results using 
the forms in Eq. (8). 

Running from Mz to Mps is done using the one-loop Standard Model RGEs given in The input values of the 
quark and lepton masses at Mz, shown in Table I, are taken from and were computed using updated Particle 
Data Group values. 





1.27±0.50 MeV 


rus/md 


19.9±0.8 


rus 


55±16 MeV 


rric 


0.619±0.084 GeV 


rub 


2.89±0.09 GcV 


mt 


171. 7±3 GeV 


rUe 


0.486 570 161±0.000 000 042 MeV 




102.718 135 9±0.000 00 92 MeV 


rrir 


1746.24±0.20 MeV 



Table I. 



The input values of the CKM angles are taken from [ll[: su = 0.2243 ± 0.0016, S23 = 0.0413 ± 0.0015, and 
si3 = 0.0037 ± 0.0005. The leptonic angles are taken fronTllil: 9soi = 33.9° ± 2.4°, daun = 45° ± 10°. The quark and 
lepton masses at Mps that result from the running are shown in Table II. 
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mass with error bar 


fractional error (%) 




0.571±0.24 MeV 


42 


nis/nid 


18.9±0.8 


4.23 


nis 


25.387±8.0 MeV 


31 


nic 


0.278±0.042 GeV 


15 


nib 


1.186±0.05 GeV 


4.2 


rrit 


86.926±4 GeV 


4.6 


nie 


0.488848231±0.000 000 042 MeV 


10-5 




103.199 06 11 ±0.000 00 92 MeV 


10-5 


nir 


1754.46±0.20 MeV 


10-4 



Table II. 



The quark mixing angles at the Pati-Salam scale are, Si2 = 0.2243 ± 0.0016, S23 = 0.0464 ± 0.0015 and S13 = 
0.0041 ± 0.0005. 

Note that wc fit m,d/m,s, which is relatively well-known from current algebra, rather than fitting m,d itself. For 
several quantities, namely the charged Icpton masses, the mass of the u quark, and the atmospheric neutrino angle, 
we will add a "theoretical error bar" to the experimental error bar in doing the fit. In the case of the charged 
leptons, we add a 1% fractional error to their masses, simply because we do not expect the forms in Eq. (8) to be 
more accurate than that. (We have made approximations of that order in deriving them.) 

In the case of m„, we add a theoretical error because the mass matrices we are using to do the fit include only 
tree- level and one- loop effects. A two- loop contribution to the 11 element of Mjj would be expected to be roughly of 
order (j^) ~ 3.5 MeV. Thus we take the prediction of the model to be that to° = ± 3.5 MeV. In other words, 
the theoretical error for m„ is about 600% of the actual value of m„, so it has no effect on the fitting. 

In the case of the atmospheric neutrino angle, there is an inherent uncertainty in the prediction of this model, due 
to the Majorana mass matrix of the right-handed neutrinos Mn being unknown and not predicted by the model. 
Because the Dirac neutrino mass matrix M^r is (to one- loop order) given by the form in Eq. (8) , which has vanishing 
first row and column, it follows that the light neutrino mass matrix, given by the usual type-I see-saw formula 
My = —MjjM^^Mpf, also has vanishing first row and column. Thus, the unitary matrix Uv required to diagonalize 
is simply a rotation by some angle Oi, in the 23 plane. From the form of Mjv, one expects generically that 
Qu = 0{e). Thus the mixing matrix of the neutrinos is given by 



/l 



MNS 











= UlUI = Ul cose^ -sin6»y 
sin cos 6v J 

where Ul is the unitary rotation of the left-handed charged leptons required to diagonalize Mj,. Consequently, 



(9) 



sin eatm = COS 9^ sin 9^^^ - sin 6^, cos O^^i cos 6l^t„ , 

= sin 6'ft„ - 0(e) ^^^^ 
sin = cos sin 9^^ — sin 9^ sin 0^^; , 
sin^so; = sin^f^;, 

where sin^^j^^, sin^^^j, and sin^fg are the 23, 12, and 13 elements of Ul respectively. In performing the fit to the 
data, we take as the neutrino mixings predicted by the model sin^^^;, sin^^^^, and sin^fg, i.e. the angles obtained 
from diagonalizing the charged lepton mass matrix. These are what are given under "model" for these quantities in 
Table III. However, we take account of the unknown 6^ by including in the "experimental error" for smOatm in Table 
III a "theoretical error" corresponding to an angle of e = 0.19 rad. 

In fitting, one has to take into account that the parameters appearing in Eq. (8) are, in general, complex. Because 
of the freedom to redefine the phases of the quark and lepton fields, most of the phases can be "rotated away" from 
the low energy theory. We will consider two cases. If the group theoretical factors denoted / and Jh are real, then 
there are two physical phases in the mass matrices of Eq. (8), which one can take to be phases of the parameters e 
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and Sh- We will denote these by 0^ and Oh- If / and Jh arc complex, then there are two additional phases, which 
we will denote 6' and 9fH- The phase 0' comes into the subleading terms of the 23 and 32 elements of and Mi^. 
The phase 6fH comes only into the 22 element of Ml- The phase 6' has only a small effect on the fit, and 0fH has 
almost none. 

Altogether, then, we fit using 9 real parameters (My, Md, Ci, C2, e, S, Sh, /, fn) and two or four phases {0^, 
6h, and if / and Jh are complex then also 0' and O/h)- With these we fit sixteen quantities: the 9 masses of the 
quarks and leptons (excluding the neutrino masses, which depend on the unknown Mr), the 3 CKM angles, the 1 
CKM phase, and the 3 neutrino mixing angles. 





model (at Mps) 


expt. (at Mps) 


off (%) 


expt. error* (%) 


rrie 


0.0004900 


0.0004888 


0.027 


1.0* 




0.1031 


0.1032 


-0.13 


1.0* 


rrir 


1.756 


1.754 


0.07 


1.0* 


ruu 





0.000571 


100 


600* 


rric 


0.342 


0.278 


23.0 


15.1 


mt 


87.24 


86.93 


0.36 


4.6 


ms/md 


18.68 


18.90 


-1.14 


4.23 


TUs 


0.0358 


0.0254 


40.8 


31.5 


TUb 


1.17 


1.186 


-1.29 


4.22 


me / rrid 


0.0886 


0.0895 


-0.99 




mb/rrir 


0.667 


0.676 


-1.35 




Vus 


0.2243 


0.2243 


0.002 


0.71 


Kb 


0.0456 


0.0463 


-1.51 


3.24 


\Vub\ 


0.00368 


0.00432 


-14.8 


11.6 


Sl3 


0.887 


0.995 


-10.8 


24.12 


sin 9soi 


0.518 


0.559 


-7.33 


7.51 


sin 0atm 


0.891 


0.707 


26.1 


28* 


sin 6*13 


0.014 


< 0.178 








7.2 









Table III. 



The results of the fit assuming / and fn are real are given in Table III. The asterisks in the "experimental error" 
column are reminders that for certain entries a "theoretical error" is included, as explained above. The masses are all 
in GeV. The CKM phase ^13 is in radians. One notes that most quantities arc fit excellently. The least good fits are 
to mc, ms, and \Vub\- Considering that 11 parameters are fitting 16 quantities, the of 7.2 is quite reasonable. 

The parameter values for this fit are e = 0.189, Ci = 1.03, C2 = -1.51, / = 0.566, fn = 0.208, 1677^^ = 2.22, 
\Qtt^5h = 2.66, 9c = 1.52 rad, 0h = 0.514 rad. Note that all these quantities are of order one. In other words, no 
small dimensionless parameters are needed to fit the quark and lepton mass hierarchies in this model. The scales 
called mu and niD in Eq. (8) are given respectively by 86.9 GeV and 0.79 GeV. The large ratio of these scales is not 
explained by the structure of the model or by symmetry, and presumably comes from the details of the sector that 
breaks the weak interactions. 



12 





model (at Mps) 


expt. (at Mps) 


off (%) 


expt. error* (%) 




0.0004888 


0.0004888 


-0.012 


1.0* 


mix 


0.1032 


0.1032 


-0.01 


1.0* 


nir 


1.755 


1.754 


0.02 


1.0* 


m„ 





0.000571 


100 


600* 


rUc 


0.317 


0.278 


14.14 


15.1 


m,i 


87.22 


86.93 


0.33 


4.6 


ms/iiid 


18.44 


18.90 


-2.39 


4.23 


rus 


0.0346 


0.0254 


36.17 


31.5 


nib 


1.17 


1.186 


-1.67 


4.22 


Trie / rrid 
' 'TIS 


0.0874 


0.0895 


-2.39 




m.h/m.r 


0.66-5 


0.676 


-1.7 




Vus 


0.224;-! 


0.224;] 


-0.018 


0.71 


Vcb 


0.0458 


0.0463 


-1.13 


3.24 


\Vub\ 


0.00382 


0.00432 


-11.5 


11.6 


(5l3 


0.889 


0.995 


-10.6 


24.12 


sin dsoi 


0.499 


0.559 


-10.7 


7.51 




0.895 


0.707 


26.7 


28* 


sin 6*13 


0.015 


< 0.178 








6.0 









Table IV. 



The results of the fit assuming / and fn are complex are given in Table IV. The parameter values for the fit in 
Table IV are e = 0.182, Ci = 0.997, C2 = -1.60, / = 0.573, = 0.224, IStt^^ = 2.16, 1&'k'^5h = 3.22, 6^ = -0.554 
rad, 9h = —1.56 rad. 

Comparison of Tables III and IV shows that the inclusion of the phases 6' and 6fH makes very little difference to 
the fits. This is not surprising, since 0' appears only on sublcading terms in the mass matrices, and OfH appears on 
the very small entry fn- For the two fits, the values of the real parameters hardly changes. The phase angles Of, 
and 6h both change by -2.07 rad, but that is essentially due to a rephasing: a shift in these two phases by a certain 
amount can be compensated by a shift in 9' together with a change in the phase of two small sublcading terms. In 
other words, what is really changing a lot between the fits in Tables III and IV is the phase 9' (which is, of course, 
zero for the fit in Table III). What this shows is that the fit is hardly affected by large changes in 9' . 

In this model there is a relation between the atmospheric angle and ^13, which is given in Eq. (10). Using the 
best-fit values given in Table III, Eq. (10) yields 

|sini9atm| = |cos6'^(0.891)-sin6i„(0.396)|, 
|sin6'i3| = |cos6'^(-O.O14)-hsin0^(O.518)|. 

If the parameter sin^j, were a real number, these equations would give a prediction for ^13 in terms of 9atm- For 
values of the atmospheric angle near maximal mixing, i.e. Oatm — 7r/4, the prediction for the 13 angle would be 
approximately given by 

I sin^igl = 0.160 - 0.72{sm9at^ - 1/V2). (12) 

However, in fact, the parameter sin^i, can be complex. Therefore, the smaller of the two values that is obtained for 
I sin 1 by solving Eq. (11) with real sm9i, is a lower bound. So, if the atmospheric mixing angle is near maximum, 
there is a lower bound on sin0i3 given by Eq. (12). 

VI. CONCLUSIONS 

The model that we have studied here is the first predictive grand unified model with a radiative fcrmion mass 
hierarchy. In a number of ways, it is as economical as a model of quark and lepton masses can be. The masses and 
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mixings of the second and third families come from only three effective Yukawa operators, shown in Eq. (1). These 
operators account for many features of the light fermion spectrum: (1) the fact that Vet is of the same order as mg/nii, 
and m^/rriT', (2) the fact that rric/mt is much smaller than those ratios; (3) the largeness of the atmospheric and 
solar neutrino angles; (4) the smallness of the 13 angle; (5) the rough equality of m° and m°; and (6) the Georgi- 
Jarlskog factor of about 1/3 between nis/mh and m^/TO,-. The masses and mixings of the first family (except for the 
solar neutrino angle) come from loop diagrams. It is remarkable that one of these loop diagrams (Fig. 1) is present 
automatically, whereas the other (Fig. 2) requires only a single additional Yukawa term to be postulated. 

It is striking that no small parameters are needed in this model to account for the dramatic hierarchies in the quark 
and lepton masses. The model yields a definite relation between the atmospheric angle and the angle ^13. 
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